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Abstract 



Q 

The possibility of realization of the phenomena of complex analytic dynamics for the re- 
alistic physical models are investigated. Observation of the Mandelbrot and Julia sets in the 
parameter and phase spaces both for the discrete maps and non-autonomous continuous sys- 
tems is carried out. For these purposes, the method, based on consideration of coupled systems, 
demonstrating period-doubling cascade is suggested. Novel mechanism of synchronization loss 
in coupled systems corresponded to the dynamical behavior intrinsic to the complex analytic 

q ■ maps is offered. 

Q\ ' 
O 

in 

p ; 1 Introduction 



It is known p!J |2], that complex analytic dynamics (CAD), studying behavior of complex maps, 
! includes a lot of interesting phenomena, for example, presence of fractal Mandelbrot and Julia sets 
^ . in the parameter and phase spaces. 

Let us start with a quadratic logistic map 

z'^X-z 2 , (1) 

where A is a complex parameter, and z is a complex variable. By definition, Mandelbrot set (fig. 1) 
is a set of points on a plane of complex parameter A, for which the orbit of an extremum z = 
of the map during iteration procedure does not escape to infinity. The Mandelbrot set contains 
the so-called "Mandelbrot cactus" (designed on figure by gray color); this is a set of points in the 
parameter plane, for which the trajectory starting of the extremum of the map converges to a periodic 
attractor. 

"Mandelbrot cactus" consists of a big cardioid, corresponding to the existence of an attracting 
fixed point, and an infinite number of "leaves", corresponding to existence of attracting cycles of 
different periods. For example, the leaves of the doubled periods are placed along a real axis. 
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Figure 1: Mandelbrot set (a) and Julia sets for the quadratic complex map with different values 
of parameter: A = 0.5 (b), A = 0.8 (c), A = 1.42 (d), A = 0.5 + 0.7i (e), A = 0.123 - 0.745i (f), 
A = 0.0315 -0.7908i (g), A = -0.282 + 0.530i (h), A = 1.16 + 0.25z (i). The gray color designates the 
regions corresponding to existence of periodic dynamics (periods are marked by respective numbers); 
the black color designates points, at which the restricted in a phase space chaotic dynamics is 
implemented; the white color means the escaping of trajectories to infinity. 

The sequences of other period m-tupling (period-multiplication) bifurcations also can be found. In 
particular, the accumulation points for period-tripling and period-quadrupling bifurcation cascades 
have been investigated in the work of Golberg, Sinai and Khanin [Hj. In the works of Cvitanovic and 
Mirheim |3J E] the universal properties of many other bifurcation cascades were studied. 

A bifurcation, which is responsible for originating a "leave", corresponds to a stability loss of a 
"parent" cycle, characterized by a complex multiplier with unit modulus and rational argument (in 
relation to 2n). If the argument of a multiplier at the stability loss is irrational, then the domains with 
fractal boundaries, filled by invariant curves arise in the phase plane, the so-called Siegel disks [H10IB!- 

Fractal pattern close to the "Mandelbrot cactus" and denoted by black color in Fig. 1, corresponds 
to existence of the chaotic dynamical regimes in the phase space. 

In figures l(b-i) the Julia sets for different values of complex parameter A are shown. The Julia 
set is a border between basins of attraction to infinity (white color) and to a periodic motion (gray 
color) in a plane of complex variable z. One can distinguish the following types of Julia sets: 

• for values A, belonging to the "Mandelbrot cactus", the Julia set is connected and enclose an 
interior basin (figs. l(b,c,f-i); 

• for values A, at which chaotic dynamics exists, the Julia set is also connected, but has no inner 
region (fig. Id); 

• for values A, outside the Mandelbrot set, the Julia set is disconnected (fig. 1/). 
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It is obvious, that ID complex map can be represented equivalently by a 2D real map (for this 
purpose it is necessary only to separate real and imaginary parts of the equation). However, the 
mentioned phenomena of CAD are intrinsic only to a very special class of the real 2D maps, namely 
for the analytic maps, obeying the Cauchy-Riemann conditions. Violation of the analyticity leads to 
drastic changes of the dynamics of the map (HJIIIIllIIllEllinj ■ Thereby, a following problem arises: Is it 
possible to specify actual physical systems demonstrating phenomena of CAD? Recently, this problem 
attracts great attention. The physical applications of complex dynamics for such problems, as the 
renormalization group approach in the theory of phase transitions and the theory of a percolation 
were discussed [HJ E3 HE1 Ej ■ In the paper of Beck [THj a theoretical possibility of the construction 
of the physical system, in which the Mandelbrot set would arise, was considered. The suggested 
approach is based on analysis of a motion of a charged particle in a double-peak potential with 
non-linear damping. The particle is driven by magnetic field, depending on time and on the particle 
velocity, and effected by external shot pulses, time intervals between which also depend on the particle 
velocity. 

In present work, we suggest a simpler and universal approach of constructing models manifesting 
the Mandelbrot set and other phenomena of CAD, which may be designed as realistic physical 
systems. This method allow us to carry out a physical experiment and present the first observation 
of the Mandelbrot set [THJ - The special structure of the Fourier spectrum of signal generated by 
experimental system at the period-tripling accumulation point is presented in 20J. 

Our method is based on using of specially symmetrically coupled identical systems, demonstrating 
transition to chaos through period-doublings. As it is known, such behavior is peculiar for a very wide 
class of nonlinear dissipative systems of various physical nature. The special kind of offered coupling 
provides a special symmetry in the multi-dimensional system, which is necessary for implementation 
of the analyticity conditions. It is a simple problem to construct the system with such coupling in 
comparison to the system, suggested by Beck. 

In section 2 the procedure of construction of the coupled systems demonstrating phenomena 
of CAD at the example of the discrete logistic maps is considered. In sections 3 and 4 we apply 
developed method to the various systems of Feigenbaum universality class, namely to the Henon 
map and to the nonlinear non-autonomous oscillator. The relation of the CAD phenomena to the 
phenomena of generalized partial synchronization is discussed in section 5. 

2 From complex quadratic map to coupled logistic maps 

Let us start with the notion that one-dimensional complex quadratic map is equivalent to the system 
of two real coupled quadratic maps with a special type of coupling. 

Separation of real and imaginary parts in the complex equation yields 



re 




irn 




(2) 



Next, we introduce the following designations 



^1 — \re + 




(3) 
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As a result we obtain a system of two coupled quadratic maps 



x[ — ► Ai — x\ + e(x 2 




(4) 



where e = (1 + /3 2 )/4/? 2 is the parameter of coupling. Note a special type of coupling in these 
equations. It can be interpreted as an identical simultaneous shift of control parameters in both 
partial systems, proportional to the squared difference of dynamic variables at each step of discrete 
time. 

It is worth nothing that the coefficient e = (l+f3 2 )/4f3 2 for any j3 is larger than 1/4. Nevertheless, 
formally we can investigate system (j3J) with any e. 

In fig. 2 we present the charts of the parameter plane (Ai, A2) for the coupled maps (jlj) at several 
values of parameter e. One can see the usual Mandelbrot set, rotated by 45°, takes place at e — 0.5. 
For 0.25 < e < +00 we have a distorted Mandelbrot set on the parameter plane. The cactus leaves of 
this Mandelbrot set correspond to existence of periodic motion of different periods. At e = 0.25, the 
set on the parameter plane, for which the point starting from the origin does not escape to infinity, 
looks like a set of strips, where the period doubling cycles occur. At s < 0.25 it transforms to a 
rhombus-like structure. At a particular e = (uncoupled logistic maps) it is a square. 

The generalization of coupling to the case of —00 < e < +00 corresponds to the original map (JTj), 
variable and parameter of which are so-called two-component numbers |21[ E2J EH EH E3 EE] • This 
is a special algebraic system, which elements are defined as follows 



According to |2T], there are three special cases: i 2 = — 1 - the usual complex numbers, i 2 = +1 
- the so-called perplex numbers, i 2 = - the dual numbers. All other algebraic number systems are 
isomorphic to complex, perplex or dual numbers depending on, whether the value of (a + b 2 )/Ab 2 
is positive, negative or zero, and are known as elliptic, hyperbolic or parabolic number system, 
respectively. In terms of parameter e these conditions look as follows: 

1) the case e > 0.25 corresponds to elliptic numbers isomorphic to complex numbers, implemented 



2) the case e < 0.25 corresponds to hyperbolic numbers isomorphic to perplex numbers, implemented 

at e = 0; 

3) the case e = 0.25 corresponds to parabolic or dual numbers. 

Thus, the existence of three topologically different structures on the plane of parameters (Ai, A2), 
namely, fractal structure, similar to Mandelbrot set, rhombus-like structure and system of strips, is 
explained by existence of three different algebraic systems of numbers. 



z = x + iy, i 



■2 



a + ib, where a, b e R. 



(5) 



at e — 0.5; 
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Figure 2: The charts of the parameter plane (Ai, A2) for the coupled logistic maps (JIJ) with different 
values of parameter of coupling: e = —0.1 (a), e — 0.0 (b), e — 0.1 (c), e — 0.25 (d), e = 0.3 (e), 
£ = 0.5 (f), e — 1.0 (g). The figures (a-c) correspond to hyperbolic numbers, figure (d) - to parabolic, 
and figures (e-g) - to elliptic numbers. 
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3 Coupled Henon maps 



The more realistic model, than the logistic map, is the Henon map j2x| |2H| 

x' — > A — x 2 — d ■ y, 

(6) 

y 

which is two-dimensional reversible map and, therefore, can be realized as the Poincare cross-section 
of some continuous system with 3-D phase space - the minimal dimensionality ensuring an possibility 
of nontrivial dynamics and chaos. 

In present section we show, that it is possible to observe the phenomena of CAD in the system of 
two coupled Henon maps. (In [29J it is shown, that phenomena of CAD such as period multiplication 
cascades survives in complex Henon map) For this purpose we carry out the following discourse, 
which starting point is the complexification of the individual map (jBJ). 

Let us consider dynamical variables x and y and driving parameter A responsible to transition to 
chaos through period-doublings as complex, that is x = x re + ix im , y = y re + iyi m , A = A re + iAj m . 
Parameter d we consider as real. Separating real and imaginary parts, we obtain: 



x re > A re x re -|- x^ m d • y T ei 

"^im * Aj m 2,X re Xi m d • yimi 



(7) 



Vim ~^ x im- 



Let us enter new variables and parameters 

X\ 2 %re P^irri} 

2/1,2 = Vre ± PVim, ^ 

Ai,2 = A re ± /3Aj TO , 
£ = (l + /3 2 )/4/3 2 . 

Then, the complexified Henon map, represented in the form of coupled real maps looks as follows 

x[ — > Ai — x\ — d ■ yi + e(x2 — x±) 2 , 

x' 2 A 2 — x\ — d ■ y 2 + e(xi - x 2 ) 2 , 
y' 2 x 2 . 

At fig. 3 the charts of the plane of parameters (Ai,A2) for various values of parameter d are 
shown. As one can see, the parameter plane contains the domains similar by the shape to the 
Mandelbrot set, distorting with increasing of \d\. At the next figure the Julia sets for system are 
represented. They also have apparent similarity to usual Julia sets. The basins of attraction of the 
fixed point (fig. 4a), cycles of period 2 (fig. 4b) and 3 (fig. 4c) corresponding to the filled- in Julia 
sets, and basin of attraction of chaotic attractor (fig. 4d) corresponding to the connected dendrit-like 
Julia set are shown. All figures in present and next sections are corresponded to the value of coupling 
parameter e — 0.5 - the value, which is equivalent to usual complex numbers and usual Mandelbrot 
set. 
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Figure 3: The Mandelbrot set for coupled Henon mappings (JHJ) with the following values of parameter 
d: —0.5 (a); —0.3 (b); 0.3 (c); 0.5 (d). Parameter of coupling e = 0.5. 
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Figure 4: The Julia sets for coupled Henon mappings with d = —0.3 for the following values of 
parameters Ai and A2: Ai = A2 = 0.1 (a); Ai = A2 = 0.5 (b); Ai = —0.8, A2 = 0.5 (c); Ai = A2 = 1.05 
(d). 
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4 Coupled nonlinear oscillators 

One of the most universal models, suitable for description of a number of real physical systems, is 
the nonlinear oscillator. Let us consider a quadratic oscillator with damping and harmonic driving 
force 

x + r yx + X + x 2 = F cos uit, (10) 

where x - dynamical variable, A - parameter of nonlinearity, 7 - parameter of damping, F and uj - 
amplitude and frequency of external driving force. It is known, that such systems demonstrate the 
opportunity of transition to chaos through a cascade of period-doubling bifurcations, for example by 
changing of parameter A with fixed F, uj and 7 [3111110] (see fig. 5a). 

Let us construct system of two coupled oscillators demonstrating phenomena of CAD, operating 
the same scheme as in the previous sections. At first, we complexify the equation of a quadratic 
oscillator in such way, that driving parameter A and variable x are complex, and then, we make 
following designations 

X\ 2 %re i f3Xi m) 
Al,2 — A re ± f3\i m , 

e = (l + {3 2 )/A(3 2 . 
As a result we obtain the system of coupled oscillators 

x\ + 7x1 + Ai + x\ — e(x2 — x\) 2 = F cos out, 
x 2 + 7x2 + A2 + x\ — e{x\ — X2) 2 = F cos out. 

At fig. 5b the chart of a plane of parameters (Ai,A2) is exhibited. Although, the represented 
picture differ by the shape from the Mandelbrot set for coupled logistic maps, but it have its basic 
properties, such as presence of leaves of every possible periods and their self-similar organization. 
The difference is determined by the existence of the region of bistability for the quadratic oscilla- 
tor (see fig. 5a). 

Analogy between basins of periodic attractors for coupled oscillators represented at fig. 6 and usual 
Julia sets is also obvious. With Ai = A 2 = —0.8 the basin of one attracting point is realized (figs. 6a). 
Then, it splits up to two basins of different fixed points (fig. 6b). With Ai = A2 = —0.575 (fig. 6c) 
one can see, that one of the basins destroys and with Ai = A2 = —0.5 there is only one attracting 
fixed point on the phase plane (fig. 6d). Figs. 6e-f demonstrate the basins of attracting cycles of 
period 2 and 3. It is necessary to note, that for non-autonomous coupled oscillators there can be a 
number of basins of attraction of various periodic and chaotic motion on the phase plane (21,22)- 
Thus, only one of them, posed at a central part of figures is associated with Julia set corresponded 
to the Mandelbrot set of figure 5b. 

Figure 7 represents the phase planes (xi,ii), (2 2 ,x 2 ) with the same values of parameters as for 
the fig. 6. 

Thus, the strategy of construction of the flow system manifesting CAD phenomena is follow (see 
scheme at fig. 8): One should take two identical elements (describable by the continuous differential 
equations), make special coupling between them and synchronize them by common periodic external 
driving force. 

8 



(12) 



■OS 

-0.8 X -0 2 -0 3 \i -02 

(a) (b) 

Figure 5: (a): Bifurcation tree for the quadratic oscillator (fTT)|) at the plane (X n , A), where X n are the 
values of dynamical variable x at the Poincare cross-section t = 2rm (n = 1,2,...). Cascade of period- 
doubling bifurcations and region of multi-stability is visible, (b): The charts of dynamical regimes on 
a plane of parameters (Ai, A2) for system of coupled oscillators (JI2J) with 7 = 0.2, F = 0.23, u = 1, 
e = 0.5. 




Figure 6: The charts of the phase space cross-section (£1,22) for coupled quadratic oscillators with 
following values of parameters: Ai = A 2 = —0.8 (a), Ai = A 2 = —0.6 (b), Ai = A 2 = —0.575 (c), 
\ 1 = A 2 = -0.5 (d), Ai = A 2 = -0.3 (e), Ai = -0.74, A 2 = -0.45 (f). The basins of attraction of the 
fixed point (figs, a-d), and cycles of period 2 (fig. e) and 3 (fig. f) at the Poincare cross-section are 
represented. 
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Figure 7: The charts of the phase plane (xi,xi) and (£2,^2) for coupled quadratic oscillators with 
following values of parameters: Ai = A2 = —0.8 (a), Ai = A2 = —0.6 (b), Ai = A2 = —0.575 (c), 
\ 1 = A 2 = -0.5 (d), Ai = A 2 = -0.3 (e), A x = -0.74, A 2 = -0.45 (f). 
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Figure 8: Scheme of construction of the physical system with CAD phenomena. 
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Figure 9: Charts of the phase plane (x 1 ,x 2 ) for the coupled logistic mappings with e = 0.5 for the 
following values of parameters Ai and A 2 : Ai = A 2 = 0.5 (a); Ai = A 2 = 0.8 (b); Ai = 0.868, A 2 = 0.622 
(c); Ai = 0.248, A 2 = 0.812 (d). The dots designate the attractors, which basins are marked by the 
gray color. The figures for a fixed point (a) and cycles of period 2 (b), 3 (c) and 4 (d) are represented. 

5 Phenomena of CAD and mechanisms of synchronization 
loss of coupled systems 

The investigation of phenomenon of synchronization loss in coupled systems demonstrating transition 
to chaos through period-doubling bifurcations is one of the mostly interesting problems of nonlinear 
dynamics and has a great fundamental and applied importance for different fields of science and 
technique. It is necessary to note, that the set of domains of periodic dynamics corresponded to 
leaves of Mandelbrot cactus, can be considered as the region of generalized partial synchronization. 
By generalized partial synchronization we mean the dynamical state of a system, in which the value 
(x\ + x\) l l 2 is bounded and varies periodically. It means, that the trajectories of one subsystem do 
not escape far from trajectories of other subsystem. At fig. 9 the charts of phase planes for coupled 
logistic maps with various values of Ai and A 2 are represented. At these pictures, the bold dots 
indicate the periodic attractors, which basins of attraction are marked by the gray color. It is easy 
to see, that in the case Ai = A 2 there is full synchronization in the system (attractor takes place at 
the diagonal line). In the case, when the parameters of subsystems Ai and A 2 are not equal, but 
belongs to the " Mandelbrot cactus" , the generalized synchronization of subsystems exists (periodic 
attractor takes place at the neighborhood of a diagonal line). 

Let us demonstrate the idea of generalized partial synchronization also with an example of phase 
portraits of coupled oscillators. One can see, that with identical values of parameters Ai and A 2 
(see fig. 8a-e) the values of variables in subsystems are coincided, that corresponds to full synchro- 
nization. In a case X 1 ^ A 2 (fig. 8f,g), when the phenomena of CAD can be realized, the values 
of dynamical variables of subsystems do not coincide, that corresponds to realization of generalized 
partial synchronization. Thus, the phenomena of CAD such as period-tripling bifurcations can be 
implemented, when the point (Ai, A 2 ) belongs to region of generalized partial synchronization. 
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6 Conclusion 



In present work we have offered a novel universal method for obtaining of the phenomena of CAD in 
the realistic models of physical systems. It is shown, that complexified system can be represented as 
two coupled real systems. Such representation is useful for some reasons. At first, often it simplify 
equations and construction of a real physical system. Secondly, with entering of special coupling to 
the system of two identical devices of any nature demonstrating the period-doublings cascade, one 
may expect the whole system to demonstrate phenomena of complex analytic dynamics. 

The coupled logistic maps, coupled invertible Henon maps and coupled nonlinear oscillators 
with periodic driving have been considered. Thus, we have studied both discrete maps, and non- 
autonomous system with continuous time. 

It has been shown, that realization of the Mandelbrot set in the parameter space of the systems 
with discreet time and non-autonomous periodically driven systems is easy problem enough. It is 
necessary to consider two elements demonstrating period-doublings with coupling, arising from the 
complexification of the variables and parameter responsible for the period doublings in the original 
system. The special interest is attracted by the problem of realization of phenomena of CAD in the 
autonomous continuous systems. 

Besides it is necessary to mark connection of the phenomena of CAD with a problem of syn- 
chronization. It is shown, that Mandelbrot set corresponds to the domain of generalized partial 
synchronization of coupled systems. Thus, the possibility of new nontrivial features of dynamics 
corresponded to synchronization of complex systems has been found. 
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